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CLIMB  PaTH  FOR  Idlv’T  ELAPSED  TIME 


The  purposes  of  this  note  are  to  examine  the  suitability  of  an  engineer¬ 
ing  method  which  lias  been  in  use  for  determining  the  path  a  jet-powered  air¬ 
plane  should  follow  to  achieve  a  specified  speed  and  altitude  in  the  least 
elapsed  time,  to  indicate  the  conditions  under  which  the  method  becomes  un¬ 
realistic,  and  to  suggest  an  alternative  procedure.  Ideaiiy,  the  optimum 
acceloration-climb  path  would  bo  obtained  as  the  solution  of  a  variation  pro¬ 
blem;  this  problem  has  been  written  out  in  very  general  form  by  Magnus  K. 
Hestenes  (RAND  EM-100).  The  solution  of  this  problem  has  not  been  found. 

To  obtain  the  method  recommended  in  this  paper  the  variation  problem  has  been 
expressed  in  the  least  general  form  which  still  takes  account  of  the  important 
mechanisms  influencing  the  time  consumed  by  the  airplane  to  accelerate  and  climb. 
The  solution  of  the  simplified  problem  is  then  approximated  by  a  systematic 
construction.  In  the  recommended  method  the  load  factor,  flight  path  angle, 
and  the  longitudinal  acceleration  are  considered  and  the  elapsed  flight  time 
is  expressed  as  the  sum  of  two  integrals,  one  corresponding  to  level  flight 
acceleration  and  the  other  to  the  climb  at  (variable)  high  speed.  The  optimum 
path  is  at  hand  when  the  sum  of  the  integruls  is  a  minimum. 

While  this  note  was  in  preparation,  the  writer  had  the  benefit  of  several 
interesting  and  fruitful  discussions  of  the  problem  with  T.  F.  Kirkwood. 

1.  First  formulation  of  the  problem. 

The  present  engineering  method  of  determining  the  optimum  acceieration- 
climb  path  is  now  to  be  considered.  The  forces  on  the  airplane  are  depicted 
on  the  following  diagram. 


n  W  dV 
D  *  —  — 

g  dt 


T  -  thrust 
D  -  drag 
W  -  weight 
L  -  lift 

g  «  acceleration  due  to  gravity) 
t  -  time 

V  *  flight  speed  (assumed  parallel  to  T) 


W 

g 


de 

dt 


Apparently  u  dV 

T  -  D  +  W  sin  9  ♦  ^  ~ 
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hence 
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(T  -  D)V 

W 


vine,  if 


Let 


(T  -  D)V 

W 


R  *  rate  of  climb  for  zero  acceleration 

o 

■  "static  rate  of  climb" 


V  sin  6 


dh 

dt 


R 


actual  rate  of  climb 


* 


and  assume  that  along  the  flight  path 

V  -  V(h) 


Then 


V  dV  dh 
+  g  dh  dt 


o/,  V  dVv 


or 


dt 


(dh  ♦  -  dV) 

_ g - 

R 

o 


Thus  the  time  to  climb,  X,  is 


(dh  ♦  -  dV) 

II  ■  —JL 


(2) 


Now  R  -  ^  where  the  thrust,  T,  ia  a  function  of  altitude,  h,  and 

velocity,  V,  and  where  the  drag,  D,  depends  not  only  on  the  altitude  and  velocity 
but  also  on  the  tangential  acceleration  and  on  the  lift  being  furnished  by  the 
wing,  which  in  turn  i3  u  function  of  weight,  angle  of  climb,  and  the  load  factor 
(i.e,,,  normal  acceleration  or  centrifugal  force).  In  the  method  of  obtaining  the 
best  climb  path  being  scrutinized,  it  is  assumed  that  Rq  is  a  function  of  altitude 

and  velocity  only,  hence  the  influence  of  load  factor  and  angle  of  climb  are 
ignored.  It  ia  further  assumed  that  the  lift  is  equal  to  the  weight  in  the 
climbing  attitude. 
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Consider  the  numerator  on  the  right  in  (2): 
i  V  dV\  .. 

1  +  I  dhj  dh  * 

On  kinematic  grounds,  it  is  possible  for  this  expression  to  vanish,  namely 
whenever  the  sum  of  the  kinetic  and  potential  energies  of  the  airplane  is  con¬ 
stant.  This  occurs  physically  when  the  component  of  the  aircraft  weight  in 
the  flight  direction  is  sustained  wholly  by  inertia.  However,  when  it  is  assumed 
that  Rq  is  a  function  of  h  and  V  only,  the  denominator  of  (2)  does  not  necessarily 

vanish  with  the  numerator  and  according  to  that  expression  the  airplane  can 
achieve  altitude  at  the  expet  se  of  speed  in  zero  or  negative  elapsed  time. 
Obviously,  the  elapsed  time  is  to  some  cxtont  underestimated  by  (2)  whenever 


dV  . 

—  ia  negative. 


Conversely,  the  eventuality  ^“>0  entails  overestination  of 


of  the  elapsed  time  by  (2)  with  Rq  ■  RQ(h,V). 

It  has  been  shown  that  the  use  of  the  formula 


-2,V2 


(dh  ♦  —  dV) 

*  (h*V) 


(2a) 


for  the  deduction  of  optimum  flight  paths  is  unrealistic  on  the  following  counts: 
The  airplane  is  restricted  to  small  climb  angles,  the  effects  of  load  factor 
(centrifugal  force)  are  not  accounted  for,  and  the  time  consumed  in  gaining 
altitude  at  the  expense  of  speed  (airplane  partially  inertia-bom)  is  underesti¬ 
mated.  Apparently  paths  determined  by  means  of  (2a)  not  be  correct. 

The  engineering  procedure  for  solving  (2a)  is  to  construct  a  chart  of  Rq 

versus  V  with  contours  of  constant  h  using  the  characteristics  of  the  specific 
airplane,  and  to  draw  on  this  chart  a  set  of  lines  joining  the  initial  and  final 
points  of  altitude  and  3peed.  The  elapsed  time  along  each  path  is  calculated 
by  numerical  integration  of  (2a).  The  trend  of  the  results  is  a  guide  to  the 
construction  of  addition.^.  path  contours,  the  process  being  continued  until 
an  optimum  path  is  found.  Unfortunately,  the  paths  so  obtained  are  usually 
such  that  the  conditions  under  which  (2a)  is  valid  and  realistic  are  violated. 

An  improved  procedure  is  needed. 

2.  Preliminary  variation-calculus  arguments. 

Equation  (2)  is  now  exploited  for  qualitative  information.  Assume  that  all 
operation  is  at  ma:ciraum  power  so  that  the  airplane  completes  a  climb  at  speeds  less 
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h(o)  -  -  initial  altitude 


h(l)  -  h2C>h1)  -  fi"*l  altitude 

V(o)  -  Vx  -  initial  speed 
V(l)  "  V2  “  final  speed. 
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.v«  form-  to  find  function.  »C 5).  v(?>  f"lflUlng 
Th.  problem  i.  then  at.t.d  in  the  fom.  to 

(Z»)  and  such  that 
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1  (h*  ♦  7v')d5 

r  [m0,v(v7 


(5) 


la  a  minimum.  '  d  V  the  following 

.  .  flrat  variation  of  X  with  respect  to  h  an 

Upon  conduction  the  first  vanav* 

aimultaneous  equations  appear 
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th 
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y  a  Ro 
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d  R. 


o  V  _ _ o 

a  v  “  g  d  h 


-  0 
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There  are 


apparently  two  a Itornativee: 


i\_ .  o 

a)  d  V  g  a  h 


or 


b)  h'  - 


and  V1 


tRo  < 


) 


in  general  (a)  !•  f-  £U^tSS^It"h2>  h  > 

follow),  hence  option  C  )  guc^  t^at  h(l)  >  h(0)  w*th 

,  (b)  there  exists  no  continuous  h(,);  th  problem  soecified 

^t°e°cept  at  a  finite  number  of  points ,  ^^Sfdoes  Liat  a  Iowan  lilt 

w'CD  and  (3)  can^be” approached  with  arbitrary  precieion. 

of  the  elapsed  time,  which  -±  by  taxing  the  pi.=~ 

differentiabiUt^condition  S  h  and  V.  Than  if 


(8) 
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ilo  V  1^0  V  0 

d  9  d  h  ^ 


one  has  that  except  at  a  finite  number  of  points 

I 


h  -  0 

v'  ■  o 


dV  dh  _  cl 1  h' 
V  ‘  dhd5  dhh 


1  *  u-r  —  m  o  or  h 

So  when  V  -  0,  either  dh 


Also:  h  - 


i  „  dV 
0,  either  V  “  0  or  dh 


arid  when  h  -  0,  either  V  -  -  -  dh 

.  titude  occurs  at  constant  ? ,  therefore:  climb  at 

"  - (io)*  8t  „  —  - 

To  decide  (2).'™??.,^  R0(»,V)  decrees.,  with 

.^1  steps.  consider  ,  ,tatlc  r.t.  o-  clisib,  increases  vdth 

tht  ei*pse‘’ time  th* rules 

- - 1 .r.l.e  et  -<-f"  conet.ent  eltlti.de 

mt  ..he  c^lss.,.  cm, stent  »ESgi.  U-e.  »2> 

►  i  rrKt.  n1  hTi  13  skfltch®^  • 

On  the  h,v  plane ,  the  fundament*.  TIL*  *  - 

A  ^  Vt  (V2  less  than  speed,  at  h2, 

I  for  maximum  static  rate  of 

i.y.  t  1  climb) 

I _ — - - - —V 

,h#hV  olane  on  which  the  followinG  equation  holds 
There  is  one  curve  on  the  h,V  plan 


a  v-  g  ^  h 


-  0. 


M.  i»  ■/  ^  Tt.V 

V  max.  By  the  nature  of  the  K^h,  ; 

^''"’’,TFnr:  iT!?tL  nimi  — 


vsm 


x 


whan  both  t.ma  vonioh,  th.refore  the  f>£“ 

rata  of  cU»b,  for  on  that  loots  both  _o  —2 

.  .  and  h  •  V_  the  corresponding  point  at  h2; 
the  point  of  this  locus  at  and  n2c,  2c 

then  if 

V  V, 

VI  1c 


*  «rrn1 erate  at  constant  altitude  to  the  point  hlc,  vic; 
it  is  advantageous  to  accelerat 

climb  along  the  locus 


to  h  V  end  then  acceler.t.  et  con.tent  Ititud.  to  V,.  The  case  V,  <  V, 
to  2c*  2c 

is  not  considered.  „  i-.'i 

It  is  to  be  emphasized  that  the  fiight  paths  paths  are, 

p“ly  de3crib 

least  general)  realistic  formulation  of  variation  problem  for 

LelSatton  A’lLab  path  of  least  elapsed  time. 

i  ^listic  'onnulatior.of  the  climb  problem  ie 
Tn  this  paragraph,  a  least  general  r  ™  for  practical  calculation 

but^in *  conjunction'1'*! th° the* f indings^ in  Paragraph  2,  it  motivates  a  suitable  approx¬ 
imation  method. 

Referring  to  equation  (1) 

W  dV  v  x  > 

T  -  D  «■  W  sin  ©  +  g  dt  * 


Assaulting  a  parabolic  airplane  polar. 


D  ♦  2~ 

•nqb  e 


D  +  D, 


where 


parasite  drag  of  airpl*116 
induced  drag 

airplane  efficiency  factor 


^  .  ’  *>  l  _.y  u*  u  *i  fi  *  m* 


%am 


P^r* *» ; 

,  \  „  J  wr  1,  v  *  J  **  -  1*. 
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wingspan 

- 


dynamic  pressure. 


The  lift  is 


W  „  d6 

L  -  W  cos  ©  ♦  g  v  dt 


th.  second  tom  being  the  centrifugal  fores.  Thus,  suh,UtuUnE  in  CD 

Vft  2  ,  „  w  *  (12) 

t(k ,v)  -  »0U.»)  *  6  ♦  ? 1  *W  -1"6  *  « 

_  .  .  _ML  is  th.  Lduesd  drag  for  M  *  the  given  speed  and 

where  u.  •  2  . 

o  ^qb  «  _  .  ,  -  nr,ri  v  only,  as  in  (12)  > 

altitude,  men  the  thrust,  T,  is  «!»« '  ” ,(  S  1°  essentially  partial  to  the 
th.  assumption  is  *“£jj to  »*•  (12)  dlmensxonl.se  by 
local  tangent  to  the  flight  the  notation 

dividing  each  term  °y  "»  ,ul  CD  S 

r  .  Ha  -  \  cD  \  pV2s  -  c  pV2  5  cL  -  -sr 

d  o  V  WD02 

D.  2  C5  2W 

.  i0  \  W  —  -  — ~  ;  C  -  — 2"  • 

</io  •  TT  •  w  ,h2.ipV2  P  v2  «t  . 


_  T(p  .VJ- 

c  W 

Then  (12)  is  ye  2  +  V 

r'.t.v)  -  ■/>/'>  *  <4>  ’v)(c0°  8  4  7  E 

.  C.pV2  .  -^§  (cos  0  .  S)2  .  Sin  8  •  | 

1  pV 


. .  v(t)  and  e(t)  along  the  flight  path. 

.  .  /n'l  relates  the  functions  p(t),  viw»  / t  \  j.  known  (by  virtue  of 

Equation  (13;  rei^,  *  courso  known  when  either  pit;  is  v  ' 

The  altitude  history,  h(^*  ^ °  (  }  e(t)  are  known,  through  the  relation 

the  standard  atmosphere  tables. 


h  (t)  -  hx  + 


V  sin  6  dt. 


W-1  U."  & 


r"’"  : » , 

j  i  .  r  u  ^  L 


The  elapsed  tine  increment  dt  associated  with  the  increments  dV,  d©  is  found  by 
separating  variables  in  (13)  regarded  as  a  dif ferential  equation;  it  is  found  that 


gdt  - 


2Q  .  I  Q  Q  '  i-  t 

— I  oos  e  d  0  +  dV  *yu  cosedVd©  *  dV2  ♦  4  —  d©2  ( t-  C-jpV2  -  sin  ©) 


2  (T-  C^V2  -  sin  © - -2  cos2  ©) 

Pv 


Thus  the  problem  is  to  determine  continuously  differentiable  functions  p(t),  V(t), 
9(t)  such  that  the  elapsed  time,  >%,  is  a  minimum,  where 


2C2  . 


•a  Cp  *2 

rr  \  •  *>  A 


gl  “ 


•^r—  0  cos  0  ♦  V  ♦  1  4  ~~77  cos  ©V©  +  (V)  ♦  4  —  &  (“T” -  C.pV  -  sin  ©) 

pv  _ *  u  v  _  o _ I _ a 


2  (r-  C.pV2  -  sin  ©  -  — |  cos2  ©) 
1  pV^ 


dt.  (14) 


The  analysis  requires  only  piecewise  continuous  differentiability  of  p,  V,  0;  but 
the  physical  situation  deletes  the  "piecewise".  Actually  only  two  of  p,  V,  9  ace 
required  for  a  complete  solution,  the  third  being  available  through  one  of  the 
following  equations 

p  -  p(h(t))  ;  h  -  h^  ♦  j  V  sin  ©  dt 


©  »  arc sin 


.  dh  , 

P  (p) 


sin  © 


The  Euler  equations  corresponding  to  (14)  are  quite  complicated.  That  arising  from 
the  variation  with  respect  to  p  is 


\2  C0  -  2C_  .  „  0  3C„ 


I  —4  cos2  e  -  c,pV2  +  -r  v  e  cos  e  -  c.  pv“  +  — f  cos4-  ©  -r  *  sin  ©  ♦ 

p  pV2  1  p2v  1  pV2 


2c  ©2  r  C 

+  —2—  -  c,pv2  (r-  2C.pV2  -  2  sin  6)  i  -4  cos2  ©  (3 T-  2C  pV2  -  3ain  ©)  -T2  ♦ 

2  1  1  u* 


f%  e*  ».(  **  i  f»,  7*  7  *  -v  5  *  ■ 

„ :  ■•■Mi. 


rt  f  :  r-  T  7  ft  i 

^  Lb  b  B  i 


_  q  _ 


*  sin  9  (2t-  C-^V2 


sin  8) 


■] 


sin  8  -  T\ 

p  v  ; 


pV 


0  cos  0  ♦ 


'2 

V* 


♦  U 


*2 
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—  (r- 


c^v2  - 


sin  0)  ■  0 


(15) 


The  aquations  corresponding  to  the  variations  of  V  and  9  are  of  comparable  complexity. 
Inasmuch  as  the  constants  C,,  C2,  and  the  dimensionless  thrust  function, 't',  are  not 
universal  but  are  characteristics  of  each  particular  airplane,  it  is  obvious  the 
-cneral  integration  of  the  three  simultaneous  equations  of  type  (15)  is  impossible, 
and  that  their  specific  integration  as  a  routine  calculation  is  out  of  the  question. 
This  observation  motivates  the  derivation  of  a  now  numerical  experimental  method  of 
the  type  described  in  paragraph  (l),  which  however  takes  into  account  the  important 
influences  of  climb  angle  and  load  factor.  And  one  has  reason  to  feel  that  such 
methods  constitute  the  only  practical  approach  to  the  problem. 


4  Derivation  of  recommendod  procedure. 

According  to  paragraph  (2;,  the  optimum  path  is  expected  to  consist  of  a  level 
portion  at  low  altitude  during  which  the  airplane  accelerates  to  a  speed  near  the 
desired  final  speed,  followed  by  a  climbing  portion  at  rather  high  3peed  (not 
necessarily  constant)  and  terminated  by  a  pull-up  if  the  desired  final  3peed  is  less 
than  the  climbing  speed,  or  by  a  levelling-out  if  the  final  speed  is  to  be  greater 
than  the  climbing  speed.  The  elapsed  tines  during  the  level  and  climbing  phases 
are  now  separately'  expressed,  and  the  optimum  path  is  at  hand  when  the  sum  of  the 
expressions  is  a  minimum. 

• 

On  the  low-and-level  phase,  It  is  assumed  that  9  -  0  *■  0,  hence,  according 

to  equation  (13)  the  time  consumed  in  level  acceleration  from  the  initial  velocity 
to  the  unknown  initial  climbing  velocity  V^c  is 


Now  the  rate  of  climb  is  V  sin  0,  hence  referring  again  to  equation  (13)  »  the 
elapsed  time  in  the  climb  phase  is 


/  ‘  dh 

1  V 

*W>  -  cT(p,V)  -  <fio(p,V)| 

(cos  9  *  rj 

f-ij 

(17) 
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The  procedure  is  as  follows:  On  a  plane  whose  coordinates  are  altitude  and 
speed,  the  initial  anti  final  points  of  the  acceleration-climb  path  are  plotted  and 
are  joined  by  path  lines  of  the  kind  shown  in  the  diagram  below. 


The  approximate  elapsed  time  formula  is  then  integrated  numerically  along  the  legs 
of  the  paths,  and  the  trend  of  the  s\im  X^  ♦  X^  is  a  guide  to  the  construction  of 

improved  path  lines.  The  process  is  continued  until  an  apparent  minimum  is  dis¬ 
covered.  The  optimum  paths  for  different  airplanes  having  the  same  type  of  pro¬ 
pulsion  system  have  of  course  a  generic  similarity  which,  when  known,  permits  the 
construction  of  a  first-guess  path  that  is  near  the  optimum  and  thereby  minimizes 
the  computation  labor  required  for  the  solution. 

Numerical  integration  by  either  the  trapezoidal  or  Simpeon  rules  is  envisaged. 
The  calculation  of  the  integrand  of  (16)  at  selected  points  on  the  level  path  line 
is  straightforward,  but  in  the  consideration  of  the  climbing  phase,  several  diffi¬ 
culties  arise: 

a.  At  the  beginning  of  the  climb  phase  the  denominator  of  (17),  which  is  equal  to 
V  sin  0,  is  zero  because  ©  is  zero  throughout  the  level  phase.  When  9  is  zero,  so 
also  is  dh;  this  indeterminacy  must  be  removed  in  order  to  start  the  numerical 
integration.  It  i3  satisf actory  to  assume  that  during  the  first  moments 

of  the  climb  the  airplane  speed  is  constant,  sin  9  -  B,  cos  9  »  1  -  r  ,  and 

2 

that  the  rate  of  climb  (and  hence  9)  is  quadratic  in  elapsed  time.  Thus,  measuring 
time  from  the  beginning  of  the  climb 


V  ein  9  -  /3t2 


Ah 


V  sin  6  d  t 


At  the  time  t 
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V  sin  9 
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hence 


t  - 


3Ah 


V  sin  0 


Also, 


8  -  at 


so 


H  .  1  e  Ah 

0  2  6  ve 


which  implies 


6  - 


2V0* 

3Ah 


where  all  quantities  in  the  last  equation  are  measured  at  the  time  t. 
Under  the  stated  assumptions,  equation  (13)  is 


io 
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Thus  at  time  t 
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2ve2 
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and  the  time  consumed  by  the  maneuver  is 
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3Ah 

V9 


(18) 


(19) 

(20) 


It  is  convenient  to  soloct  a  Ah  of  about  500  feet;  all  other  quantities  in  these 
expressions  are  then  known  in  terms  of  the  known  V  and  h.  One  is  now  in  a 
position  to  continue  the  integration  of  (17)  along  the  selected  climb  path. 

b.  At  subsequent  points  (labelled  5*  6,  7,  . ..,  j,  ...  on  the  accompanying 
figure)  along  the  path  line  on  the  h,  V  plane,  only  h,  V,  and  dh  are  definitely 

known,  but  0,  0,  and  V  are  required  as  well  in  order  to  calculS^e  the  integrand 
in  (17)  at  these  points.  Following  a  suggestion. by  Kirkwood,  these  quantities 
may  be  estimated  by  a  continuing  integration  of  0  as  follows: 
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At  point  J,  h.  V,  and  §  .re  kno-n,  and  .t  point  1  -  1.  *U  ''”»"Uti<!9 


are  known . 


Then  at  j 


I  «  i  (a  n 

dV  *1 


^  ^  • 
h  a  Valn8  -  “  V 


gativo  quantity  in  the  figure  as  drawn.) 


=>  y  .  i|  V  sin  0 


r  -  <fr 


*  <flo  («»  6  .  *  .in  6  ♦  £ 


bo  at  j  »  2  /  .u\ 

r  -  *  •rlo  (cos  6  *  f)  *  5ln  6  (L  *  e") 


e3  -  (r^~  A.  ^  ~  V 

separate  uT^/a  f£Ll 

may  be  assumed.  Then 


9,  ♦  ©i 

9,  “  + 


(V  sin  e)  ^  ♦  (V  3in  0) 


hence 


(e  -  e^v  sin  e)y+  (v  sln  e^i-i] 


The  equation  obtained  by  comparing  expressions 


(22)  and  (23)  may  be  solved  for  0^ 


by  cut-and-try,  whereupon  either  expression  may  oe 
by  (21) 

V  -  n  (V  sin  0) . 

i  J 


may  be  used  to  calculate  9y  Finally', 
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.  „  <„  replaced  bv  zero  in  the  foregoing  equations, 

When  climbing  at  constant  velocity,  n  is  replaced  Dy 

paration  and  will  constitute  a  subsequent  report. 
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